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We present an analytical computation of the asymptotic temporal behavior of the information 
geometric complexity (IGC) of finite-dimensional Gaussian statistical manifolds in the presence 
of microcorrelations (correlations between microvariables). We observe a power law decay of the 
IGC at a rate determined by the correlation coefficient. It is found that microcorrelations lead to 
^-^, the emergence of an asymptotic information geometric compression of the statistical macrostates 

C^ ' explored by the system at a faster rate than that observed in absence of microcorrelations. This 

finding uncovers an important connection between (micro)-correlations and (macro)-complexity in 
Gaussian statistical dynamical systems. 
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7^' I. INTRODUCTION 

Cd The study of complexity [Ij has created a new set of ideas on how very simple systems may give rise to very complex 

behaviors. In many cases, the "laws of complexity" have been found to hold universally, independent of the details of 
the system's constituents. Chaotic behavior is a particular case of complex behavior and it will be the object of the 
present work. In this article we make use of the so-called Entropic Dynamics (ED) [2] and Information Geometrody- 
"^^ " namical Approach to Chaos (IGAC) [3,|J1- ED arises from the combination of inductive inference (Maximum Entropy 
JTJ ■ Methods, 5, 6]) and Information Geometry [7|. ED is a theoretical framework whose objective - among others - is to 
Qs ' derive dynamics from purely entropic arguments. The applicability of ED has been extended to temporally-complex 
^sj 1 (chaotic) dynamical systems on curved statistical manifolds M.s resulting in the information geometrodynamical ap- 
PsJ proach to chaos (IGAC) 3]. IGAC arises as a theoretical framework to study chaos in informational geodesic flows 

1*-! describing physical, biological or chemical systems. A geodesic on a curved statistical manifold M.s represents the 
f**^ ■ maximum probability path a complex dynamical system explores in its evolution between initial and final macrostates. 
^^ \ Each point of the geodesic is parametrized by the macroscopic dynamical variables {6} defining the macrostate of the 
'~H . system. Furthermore, each macrostate is in a one-to-one correspondence with the probability distribution {p(X|0)} 
representing the maximally probable description of the system being considered. The set of macrostates forms the 
parameter space I?e while the set of probability distributions forms the statistical manifold M^s- IGAC is the in- 
y^ ' formation geometric analogue of conventional geometrodynamical approaches [8|, |9| where the classical configuration 
space Ye is being replaced by a statistical manifold M^s- This procedure affords the possibility of considering chaotic 
dynamics arising from non conformally flat metrics (the Jacobi metric is always conformally flat, instead). It is an 
information geometric extension of the Jacobi geometrodynamics (the geometrization of a Hamiltonian system by 
transforming it to a geodesic flow [lO|). The reformulation of dynamics in terms of a geodesic problem allows the 
application of a wide range of well-known geometrical techniques in the investigation of the solution space and prop- 
erties of the equation of motion. The power of the Jacobi reformulation is that all of the dynamical information is 
collected into a single geometric object in which all the available manifest symmetries are retained- the manifold on 
which geodesic flow is induced. For example, integrability of the system is connected with existence of Killing vectors 
and tensors on this manifold. The sensitive dependence of trajectories on initial conditions, which is a key ingredient 
of chaos, can be investigated from the equation of geodesic deviation. In the Riemannian Q and Finslerian Q (a 
Finsler metric is obtained from a Riemannian metric by relaxing the requirement that the metric be quadratic on 
each tangent space) geometrodynamical approach to chaos in classical Hamiltonian systems, an active field of research 
concerns the possibility of finding a rigorous relation among the sectional curvature, the Lyapunov exponents, and 
the Kolmogorov-Sinai dynamical entropy (i. e. the sum of positive Lyapunov exponents) [ll| . 

Understanding the relationship between microscopic dynamics and experimentally observable macroscopic dynamics 
is a fundamental issue in physics |12h15| | . An interesting manifestation of such a relationship appears in the study of 
the effects of microscopic external noise (noise imposed on the microscopic variables of the system) on the observed 



collective motion (macroscopic variables) of a globally coupled map [T^| . These effects are quantified in terms of the 
complexity of the collective motion. Furthermore, it turns out that noise at a microscopic level reduces the complexity 
of the macroscopic motion, which in turn is characterized by the number of effective degrees of freedom of the system. 

In this article, using statistical inference and information geometric techniques, we investigate the macroscopic 
behavior of complex systems in terms of the underlying statistical structure of its microscopic degrees of freedom in 
the presence of correlations. We compute the asymptotic temporal behavior of the information geometric complexity 
of the maximum probability trajectories on finite-dimensional Gaussian statistical manifolds in the presence of micro- 
correlations. We observe a power law decay of the IGC at a rate determined by the correlation coefficient. The ratio 
between the IGC in the presence and in the absence of microcorrelations is explicitly computed. We conclude that 
microcorrelations lead to the emergence of an asymptotic information geometric compression of the explored statisti- 
cal macrostates (on the configuration manifold of the model in its evolution between the initial and final macrostates) 
that is faster than that observed in absence of microcorrelations. 

The layout of the article is as follows. In Section II, we briefly discuss Gaussian statistical models in absence 
and presence of microcorrelations. In Section III, we introduce the Gaussian statistical model being considered. 
We compute the Ricci scalar curvature and the geodesic trajectories of the system. In Section IV, we compute the 
asymptotic temporal behavior of the dynamical IGC of the model. Our conclusions are presented in Section V. 

II. ON GAUSSIAN STATISTICAL MODELS AND MICROCORRELATIONS 

In this Section, we introduce the notion of Gaussian statistical models (manifolds) in the presence of correlations 
between the microscopic degrees of freedom (microvariables) of the system (microcorrelations). 

A. Statistical Models in Absence of Microcorrelations 

Consider a Gaussian statistical model whose microstates span a n-dimensional space labelled by the variables 
{X} — {xi, X2,..., Xn} with Xj e R, Vj = 1,..., n. We assume the only testable information pertaining to the 
quantities Xj consists of the expectation values {xj) and the variance A.Xj. The set of these expected values define 
the 2n-dimensional space of macrostates of the system. A measure of distinguishability among the macrostates of 
the Gaussian model is achieved by assigning a probability distribution P {X\Q) to each 2n-dimensional macrostate 

Q = |((i)0j,(2)0j)}^^ _^.^_^ = {{{xj) , Axj)}^ jj.g. The process of assigning a probability distribution to each state 
endows Ms with a metric structure. Specifically, the Fisher- Rao information metric g^^ (O) 7] is a measure of 
distinguishability among macrostates on the statistical manifold A^s, 

g^AQ) ^ J dxp {x\e)d^\og p {x\e)d,iog p {x\e) = 4 J dxd^^p {x\e)d,^p (x\e), (i) 

with n, u — 1,..., 2n and 9^ — g^- It assigns an information geometry to the space of states. The information metric 
Qfii^ (0) is a symmetric, positive definite Riemannian metric. For the sake of completeness and in view of its potential 
relevance in the study of correlations, we point out that the Fisher-Rao metric satisfies the following two properties: 
1) invariance under (invertible) transformations of microvariables {x} G X; 2) covariance under reparametrization of 
the statistical macrospace {9} G Vg. The invariance of g^^i, (9) under reparametrization of the microspace X implies 

0, 



X CW 3 X ^-^ y '^^^ f (x) e y CW =^ p {x\9) K-> p' {y\9) 



dx 



(2) 



The covariance under reparametrization of the parameter space "Dg (homeomorphic to A4s) implies [7[ 

89° d9^ 



Vg3 9^9' "=^' / (9) e Ve' =^ g^. (9) ^ g'^, {9') - 
where 



89'^" 89'" 



9afi (d) 



(3) 
e=/-i(e') 



<. {0') = / dxp' {x\9') a; logp' {x\9') 8, \ogp' {x\9') , (4) 



with d' = -T^Tjr and p' (x\9'^ = p(x\9 — f~^ [6'^^. Our 2n-diniensional Gaussian statistical model represents a 
macroscopic (probabilistic) description of a microscopic n-dimensional (microscopic) physical system evolving over 
a n-dimensional (micro) space. The variables {X} = {xi, a;2,.--, a^n} label the n-dimensional space of microstates of 
the system. We assume that all information relevant to the dynamical evolution of the system is contained in the 
probability distributions. For this reason, no other information is required. Each macrostate may be thought of as a 
point of a 2?i-dimensional statistical manifold with coordinates given by the numerical values of the expectations ^^'9j 
and ^"^'Oj. The available relevant information can be written in the form of the following 2n information constraint 
equations, 
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(5) 



The probability distributions Pj in ([5)) are constrained by the conditions of norniahzation, 
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(6) 



Information theory identifies the Gaussian distribution as the maximum ent ropy d istribution if only the expectation 
value and the variance are known [17]. Maximum relative Entropy methods [a,ll3ll9i allow us to associate a probability 
distribution P {X\Q) to each point in the space of states Q. The distribution that best reflects the information 
contained in the prior distribution m{X) updated by the information {{xj) , Axj) is obtained by maximizing the 
relative entropy, 



^(9) 



d^xPix\e)iog 



p 



m{x) y ' 



(7) 



where m (X) is the prior probability distribution. As a working hypothesis, the prior m {X) is set to be uniform since 
we assume the lack of prior available information about the system [20|. We assume uncoupled constraints among 
microvariables Xj. In other words, we assume that information about correlations between the microvariables need 
not to be tracked. Therefore, upon maximizing ([7]) given the constraints ([5]) and (|6]), we obtain 
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p{x\e) = l[p,{^x,\('^e„('^e, 



(8) 



i=i 



where 



P. 



{x,\^%S''e,) 



{2t:(j]) ^exp 



2a] 



(9) 



(2)0, "^^ Ax., = a. 



and, in standard notation for Gaussians, ^^'6j = (xj) = //.-, ^'^^Oj =' Ax^ = CTj. The probability distribution 
encodes the available information concerning the system. The statistical manifold A4s associated to ([8]) is formally 
defined as follows, 
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Xjl^i. 



(10) 



unless specified otherwise) I^ 
from (EI) is im 



and O belongs to the 2T7,-dimensional parameter space Vq ~ [I^ x X^] . The parameter space Vq 

id 
(0, + oo)^. The line element ds^ = g^^ (6) dQ^dQ" arising 



where X e 

(homeomorphic to Ms) is the direct product of the parameter subspaces T^ and X^, where (in the Gaussian case 



-oo. 



oo) and T^ 




(11) 



with /i, J/ = 1, 



2n. 



B. Gaussian Statistical Models in Presence of Microcorrelations 

Coupled constraints would lead to a " generalized" product rule and to a metric tensor with non-trivial off-diagonal 
elements (covariance terms). In presence of correlated degrees of freedom {xj}, the "generalized" product rule 
becomes, 



(12) 



Ptot {Xl,..., Xrr) = Y[ P] {Xj) ™"£4.'°"" p^^^ {Xi,..., X^i) ^ J| Pj (Xj) , 

where, 

Ptot (Xl,..., X„) = Pn{Xn\xi,..., Xn~l) Pn-l{Xn-l\xi,..., Xn-2) ■■■P2ix2\xi) Pl{xi) . (13) 

Correlations among the degrees of freedom may be introduced in terms of the following information-constraints. 



^j = fj (a^iv, a^j-i), Vj 



(14) 



In such a case, we obtain 

Ptot (Xl,..., Xn) = 6{Xn - fn {xi,..., Xn-l)) 5 {Xn-l ~ fn-1 {xi,..., X„_2)) ...5{x2 - h {xi)) Pi {xi) , (15) 

where the j-th probability distribution P, [xj) is given by. 



Pj i^j) ^ I ■■ I dxi...dxj^idxj+i...dxnPtot (a^i,---, Xn) 



(16) 



A formal manner in which correlations are introduced in probability theory is as follows. Given two arbitrary randomly 
distributed variables Xi and a;2, consider the problem of finding a linear expression of the form ci -I- C2X2, involving 
real constants Ci and £2 such that ci -1-02X2 is the best "mean square approximation" to xi. The best approximation 
is such that 



[xi - ci - C2X2) ) = min ( {xi - c\- C2X2) 

I Cl, C2 \ 

where the minimum is taken with respect to all real constants c\ and C2. To solve this problem, let 

A^i = {xi) , o\ = {{xi ~ {xx)f^ , ii2 = {^2) , <jI = ((a;2 - {X2))j 



and introduce the quantity [22], 



dof ((a;i - (a;i)) {x2 ~ {X2))) (a;ia;2) - /x^/ij 



(17) 
(18) 
(19) 



(T1CT2 0'icr2 

The quantity r is the so-called correlation coefficient of the random variables xiand X2- For the sake of convenience, 
we may introduce the "normalized" random variables. 



dcf Xi ~ fl, dof X2 — fi2 

Vi = and, 772 ^ 



The problem in (J17p can now be reduced to. 



min ((ryi - ci -02772)' 

Cl, C2 



(Jl 
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Cl, C2 



<J2 



{l^r')+cl + {r-C2Y 



= 1 - r^ > 0. 



(20) 



(21) 



The minimum is achieved for ci — and 02 — r, where r lies in the interval — 1 < r < -|-1. 

In our work, correlations among the microscopic degrees of freedom of the system {xj} (microcorrelations) are 

conventionally introduced by means of the correlation coefficients r, 



(micro) 



, N dcf \XiXj) \Xi) \Xj) 

r [Xi , Xj) = ^ ^^ — u^^LL ^ Tffith a,. 



(JiCTi 



{xi - {xi)y 



(22) 



with 



(micro) 

r. . 



€ (—1, 1) and i, j = 1,..., n. For the 2?T,-dimensional Gaussian statistical model in presence of microcor- 



relations, the system is described by the following probability distribution P{X\Q), 



p(x|e) 



[(27r)"detC(e)]- 



— exp 



(X ~ Mf ■ C'^ (6) • {X - M) 



^Y[{2i,a^) ^exp 



j=i 



2a| 



, (23) 



where X = (xi,..., a;„), M = (/^j^,..., /i„) and C(0) is the {2n x 2ri)-dimensional (non-singular) covariance matrix. 



III. THE MODEL 



In this Section we focus on microcorrelated Gaussian statistical models with 2n — 4. For n = 2, (|23l) leads to the 
probability distribution P (x, y\^J■x^ o'x, /^t,: fy) which takes the form, 



P{x, 2/|m^, <Jx, fJ.y, (Jy 



exPi-2(T^7^ 



r-j ^' ^^ 1 r-j 



2TraxCryy^l - 



(24) 



where era; > 0, a-y > 0, r £ (—1, + 1). Substituting (l24l) in ([T]), the Fisher-Rao information metric 
5mi' (Ma;> cTx, fJ-y, cTy-, r) bccomes, 



/ a^(r--l) " ^ T,<T„(r--l) 



5^1^ (A*a;> CTa:, /Xy, CTj,; r) 











\ 







J(r^-l) 







a^<7y(r'-'~l) 



crxO-„(r^-l) 



(25) 



The infinitesimal line element ds]^^ relative to g^iy (/i^., a^, fJ^y, Cj,; r) is given by, 



ds 



Ms 



511 (o-^;; ?■) d^^ + 333 (c^y; ?") ^My + 522 {(T.^; r) dal + ga (^y; r) dal + 2gi3 (cr^, cry; r) dfi^d^iy 



+2^24 (cx, CTy; r)dcrxd(J 



v 



(26) 



where, 



511 (^.; -) = -^2(,Li)' 513 (-., ^.; = ^^^^('2_i)' 522 (a.; r) = -^^^y 



''' ^^^' ""'' '^ ^ a.aj(r^-iy ''' ^^^ ''^' "^ = ^..,,('^-1)' ''' (^^' '^^ ^ "^fF^' 



542 (ctx, ct^; J') = 



'a:"j^ (?' ^ 1) 
cr:rcra (r^ - 1) 



94ii<Jy;r) = - 



r 
xuy (r^ — 1) 

2-r2 



a2(r2-l)- 



(27) 



The analytical study of the IGAC arising on a curved statistical manifold with infinitesimal line element given by 
ds^j^ in (|26p turns out to be rather difficult. Hence, as working hypothesis, we are going to assume two correlated 
Gaussian-distributed microvariables characterized by the same variance, that is we assume cr^; = Uy = a. Thus, the 
simplified line element becomes. 



dsMs = 511 (c^x; r) d^il + 533 (cr^; r) d^l + 2gi3 (cr; r) dn^d^iy + [522 (ct; r) + 344 (cr; r) + 2.g24 (c; t')] da'^ . (28) 

The new Fisher-Rao matrix g^j^^ (/i^, /i , cr; r) associated with line element ds^ in (|28p becomes, 



1 r n 

I I 7^^ 2(r^-l) " 

5m'' (A^a;J A'yi '''''') ^ ^2 I 2(r^^ ^T^~^ ^ 

4 



(29) 



We will study the information dynamics on curved statistical manifolds j^y^°"'^ '^ ^°'^'^^) g^^^j ^^o-corrc a lonsj ^j^j^ 

. f, . , , IT , / J 9 \ correlations 1/79 \ no-corrclations . , ^m t i i / 7 9 \ correlations 

mfimtesimai Ime elements (ds^gj and (afs^gj , respectively. 1 he line element (as^^j 

is defined by, 



/ 7 9 X correlations dct 1/ 1 ,2, 1 72 ^T" ,/72 



(30) 



, ., / T 9 \ no-corrclations . , . , p / ? 9 \ correlations . ,1 i. . ,1 , i 

while ^as^ 1 is obtained rroni ydsj^ \ in tne limit tliat r approacnes zero. 



A. Information Geometry of The Model 

Consider the information dynamics of the Model introduced in Section II. The Fisher-Rao line element 
{ds'j^ ) of such statistical model M.g'^° ^ '""'^ is given in (1501) . The inverse metric tensor 5'^" (/i^, fi , a; r) 

is given by, 



5'"'(Mx> fJ-y^ <T;r)=a 



( 4p-l) 2r(r^-l) \ 

r'-'—i r'^—i 

2r(r^-l) 4(r^-l) 

r^— 4 r2_4 t) 

V iy 



(31) 



The metric tensor g^ii^ (/x^, /i , cr; r) and its inverse g'^" (/i^, /i , ct; r) are necessary to determine the Christoffel 
connection coefficients T^a of the manifold J\/ih°^^°'^ lonsj^ Recall that the connection coefficients Ff. are defined as 

m. 

In our case, the non-vanishing connection coefficients are given by, 

III, rli -, lo 111 



T^3 T^3 T^3 T^l T-'l T^ 

11 ~ A „-? 1 _ ' ^ 12 ~ ^ 21 — o /„9 1 \ _ ' ^ 13 — 1^ 31 — ^ _ ' J^ 



3 



-^23 — ^32 = 1 1^33 = • (33) 

Once the non- vanishing components of F*^- are obtained, we compute the Ricci curvature tensor TZij defined as [23|, 

TZij = SfcFj^ - 9jF^j, -I- Fj^F'^'„ - F^F^-^. (34) 

Substituting p3p in p4p , we obtain the non- vanishing Ricci curvature tensor components TZij , 

2 (r^^^' '' ~ '' " "4(r2-l)^' '" " 2 (r2 - 1) a 
Finally, we compute Ricci scalar curvature TZm, {t) 



-^11 - 77712 TT:;^' -^12 - -R2I - -;j7-2 7T-2. ^22 - :r7-2 7T-2 ' ^33 - - -3 ■ (35) 



nMAr)^^'n,,g'^. (36) 



Substituting (1551) and (I5T|) in (I55|) . T?.^^, (r) becomes, 

T^A^, (?•) = .9"^ii + 2<?i2^i2 + 5''fi22 + 5''i?33 - -| (37) 

Therefore, we conclude that J\4s^^° ^ ^°^^ is a curved statistical manifold of constant negative curvature. 

B. Information Dynamics on Ms 

The information dynamics can be derived from a standard principle of least action of Jacobi type [2!] . The geodesic 
equations for the macrovariables of the Gaussian ED model are given by nonlinear second order coupled ordinary 
differential equations, 

d2e^ „„ dQ" dQP 



dr^ '^-^^-O- (^«) 

The geodesic equations in psp describe a reversible dynamics whose solution is the trajectory between an initial 
Q(initiai) g^j^jj ^ final macrostate Sf^"'^^. The trajectory can be equally well traversed in both directions. In the case 



under consideration, substituting ((29|) in ((38|) . the three geodesic equations become, 
d^fi^ (t) 2 d/z^ (r) do- (r) 
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dr^ (T (t) dr dr 

dVy (t) 2 d^j^ (r) dcr (r) 



dr^ 



a (r) dr dr 



Sa{T) 1 (da[T)V 1 1 1 fd^i^{T)Y 1 1 1 (diiyir) 



dr^ (7 (r) \ dr y 4 r^ — 1 cr (r) \ dr 

r 1 d^^ (r) d^iy (r) 

4 (r^ - 1) cr dr dr ' 



4 r2 - 1 cr (r) \ dr 



(39) 



Integration of the above coupled system of differential equations is non trivial. A detailed derivation of the geodesic 
paths is given in the Appendix. After integration of ([39]), the geodesic trajectories become, 



fJ-x ('t; ^) = - 



2(70^1 



a (r; r) — 2a, 



VMr) 1 + cxp haoy/MTjA 
(^aoy^A{ry. 
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2aoA2 



V^lrj 1 + exp hao^/M^' 



exp I 



+ exp(2aoy/MrJTy 



(40) 



where, 



dcf Al + Al-rA,A2 
^^'^ 4(1 -r2) ■ 



(41) 



Notice that for any real value of Ai and A2, < (Ai - ^2)^ == AI + AI-2A1A2 < Al + Al-rAiA2 and 4 (l - r^) > 
for r E (—1, 1). It then follows that, A{r) > 0. Note that cr(r; r) £ (0, +00) while /i^, (r; r) and Hy{T; r) £ 

(—00, +00). 

IV. THE INFORMATION GEOMETRIC COMPLEXITY AND MICROCORRELATIONS 



We recall that a suitable indicator of temporal complexity within the IGAC framework is provided by the information 
geometric entropy (IGE) Sms (t) [SjIJI' 



Sms (t) == log ^o^ 



V, 



(geodesic) 

e 



(r) 



(42) 



The information geometric complexity (IGC) is defined as the average dynamical statistical volume vol 
given by, 
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= lini — / dr^ vol 



r^-oo \ T 
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ir') 



where the "tilde" symbol denotes the operation of temporal average. The volume vol 



V, 



(geodesic) 

e 



(r') 



vol 



^(geodesic) / 






(43) 
is given by, 

(44) 



where N is the dimensionality of the statistical manifold M^ and PiMs g) (^^vi ^ ) is the so-called Fisher density 
and equals the square root of the determinant of the metric tensor g^^ (0), 



P,M.^,,{e\...,e-)'='^g{9\...,9-). 



(45) 



The integration space X>q ° ^ (r') in (J44l) is defined as follows, 



,j-»(gcodcsic^ 



(r') ='{e = 



(.'„ 



]7V 
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f{r')}, 



(46) 



where A: = 1,.., N and 6* 



^(geodesic) 



is the volume 



ei'' (s) with < s < r' such that 61* (s) satisfies (|38l)- The integration space p^^"""**"^^ (,-') 
in (j46|) is an A^-dimensional subspace of the whole (permitted) parameter space Vq . The elements of I?q ° '^^^'^ (r') 
are the A^-dimensional macrovariables {Q} whose components 9 are bounded by specified limits of integration 6 (0) 
and 6 ' (r') with k — 1,.., N. The limits of integration are obtained via integration of the geodesic equations. 
Formally, the IGE Sm^ (t) is defined in terms of a averaged parametric {N + l)-fold integral (r is the parameter) 

over the multidimensional geodesic paths connecting 9 (0) to 9 (r). The quantity vol Vq° ^^^''' 

of the effective parameter space explored by the system at time r'. The temporal average has been introduced in 
order to average out the possibly very complex fine details of the entropic dynamical description of the system on Ms- 
Thus, we provide a coarse-grained-like (or randomized-like) inferential description of the system's chaotic dynamics. 
The long-term asymptotic temporal behavior is adopted in order to properly characterize dynamical indicators of 
chaoticity (for instance, Lyapunov exponents and the Kolmogorov-Sinai dynamical entropy) eliminating the effects 

of transient effects which enter the computation of the expected value of vol Vq" 

one observes that typical initial conditions behave in an apparently chaotic manner for a possibly long time, but then 
asymptotically approach a nonchaotic attractor in a rapid fashion. 

(geodesic) 
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In the case under consideration, vol 
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(r) is given by. 
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(47) 
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) 1 J /i / \ dof / 
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Since vol 



V, 



(geodesic) 



e 



(r) must be positive by construction, it must be the case that ^1^12 < 0. Without loss of 

generality, we assume Ai — —A2 = a G M.. As a side remark, we point out that from (PU]) it follows that an increase 
(decrease) in |/i^ (r)| is followed by an increase (decrease) in |/ij,(r)|. Therefore, internal consistency (consistency 
with the class of geodesic paths considered) requires that we limit our analysis to positively correlated microvariables. 



that is r G (0, 1). Finally, the average dynamical statistical volume vol 



V\ 



(geodesic) 

e 



(r) 



becomes, 



vol 



V, 



(geodesic) 

e 



(r) 



vol 



.*(gcodcsic) 



251,81-ouesie; / /-, 



dr' 



9' 



2cto A^ (r) 



1 - exp I -(Tn^/AirJT 



(49) 



In the long-time limit, the asymptotic behavior of the IGC becomes, 

a^ g'2 (r) 1 



vol 



V. 



(geodesic) 

e 



(r) 



(50) 



cor. 2cto y^5 (r) T 

Thus, comparing the asymptotic expressions of the IGCs in the presence and absence of microcorrelations, we obtain. 



vol 



V. 



(geodesic) 

e 



(r; r) 



vol 



V. 



(geodesic) 

e 



(t;0) 



g'^{r)A^{Q) 1 /4(4-r2) / 2 + 



Ai{r)g'h{{)) 2t Y (2-2r2)H4(l-r-2) 



= y'Ms (r) 



(51) 



Written alternatively, 



..(geodesic) 



^,^^eo.e=.c; ^^. ^^ 



4(4-r2) 



2i V (2-2r2)^ \A{l~r^) 



vol 



V. 



(geodesic) 

e 



(t;0) 



(52) 



We emphasize that Fms {t) is a monotonically decreasing function of r, that is Fms (''i) — ^Ms ('^2) for any ri < r2 
with ri, r2 & (0, 1) and < J^ms (^) < 1 for r € (0, 1). We observe an asymptotic power law decay of the IGC 
in ([50|l at a rate determined by the correlation coefficient r. The ratio between the IGC in the presence and in 
the absence of microcorrelations in (j5ip leads to conclude that microcorrelations cause an asymptotic information 
geometric compression of the explored statistical macrostates at a faster rate than the that observed in absence of 
microcorrelations. Our finding presented in (|52p shows an important connection between (micro)-correlations and 
(macro)-complexity in Gaussian statistical models. 

V. FINAL REMARKS 

In this article, we presented an analytical computation of the asymptotic temporal behavior of the IGC for a finite- 
dimensional microcorrelated Gaussian statistical model. The ratio between the IGC in the presence and in absence 
of microcorrelations was explicitly computed. We observed a power law decay of the IGC at a rate determined by 
the correlation coefficient. Specifically, the presence of microcorrelations lead to the emergence of an asymptotic 
information geometric compression of the statistical macrostates explored by the system at a faster rate than that 
observed in absence of microcorrelations. This result constitute an important and explicit connection between (micro)- 
correlations and (macro)-complexity in statistical dynamical systems. The relevance of our finding is twofold: first, it 
provides a neat description of the effect of information encoded in microscopic variables on experimentally observable 
quantities defined in terms of dynamical macroscopic variables [27|; second, it clearly shows the change in behavior of 
the macroscopic complexity of a statistical model caused by the existence of correlations at the underlying microscopic 
level. 

We are confident that this work constitutes an important preliminary step towards the computation of the asymp- 
totic behavior of the dynamical complexity of microscopically correlated multidimensional Gaussian statistical models 
and other models of relevance in more realistic physical systems. In principle, our approach extends its application 
to arbitrary statistical models that may arise upon maximization of the logarithmic relative entropy subject to the 
selected relevant information constraints. In particular, our findings here presented could find practical applications 
in the statistical analysis of biological and social systems since Gaussian statistical models are of primary importance 
in statistical studies [25-|. However, our ultimate hope is to extend this approach in the field of Quantum Information 
to better understand the connection between quantum correlations (entanglement) and quantum complexity [26-30i [. 
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Appendix A: Integration of the Geodesic Equations 

From the first and second Equations in (|39p . we obtain, 



^ = 2iM.„diiW=2^, (Al) 



respectively. From (|A1[) it follows that, 

A, (r) = Aia^ (r) and /i, (r) - ^2^^ (r) , (A2) 

where Ai and A2 are real constants. Substituting (|A2I) in the third Equation of l\'S9\i we obtain, 

a (r) a (r) - d^ (r) + ^^^f^^<^' (r) = 0. (A3) 
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Therefore, the coupled system of differential equations reduces to, 

A, (r) - Aia2 (r) = 0, 

Aj, (r) - ^2^2 (r) = 0, 

a{T)a{T)-a^{T)+Air)a^{T) = 0, (A4) 

where we recall that. 

We now proceed as follows: integrate the nonlinear differential equation a (r) a (r) — &^ (r) + A (r) ct"' (t) = and 
then calculate /i^. (r) and fiy (t). 

Letting y (r) = a (r), the first nonlinear differential equation to integrate becomes, 

yiT)y{T)-f{T)+Air)yUr)=0. (A6) 

Performing the following change of variables, 

yir)^^^i{r) (A7) 

equation (jA6l) becomes 

ix -i^ +A (r) x^ = 0. (A8) 

Equation (jASp can be integrated as follows. Performing the following additional change of variables, 

dx (t) 

leads to 



= z (x) (A9) 



X = zz 



' and, X = {z"z + z'^) z, (AlO) 



with z' ^ ^. Substituting ([M]) and (|JTO)) into (IM)) . we find 

z" + Air)z^O. (All) 

Integration of (|Alip yields 

z{x) :=: A3 cos (^/MJjx^ + At sin (^/A{r)x^ , (A12) 

where ^3 and A4 are reaZ constants. Recalling that x = ^^ = z (x), we have 

'^'^ = fdT + A^. (A13) 



A3 cos (^/A(r)x^ + A4 sin (^/A(r)x] 
Reality conditions imply that A3 — and without loss of generality we can set A^ = 0. Integration of (jA13|) leads to. 



1 l-exp[2Ai^Air)T] 

X (r) = — ^^ arccos ) ( I . (A14) 

V -^ (r) \l + exp {2A4 y/AJZ 



Finally, recalling that y (r) — dP' ^^"^ 2/ (''') ^ "' (''')' '^^ S^* 



exp { Ai^y A ir)T 

o-(t)=2A4 ^7 ^. (A15) 

l + expf2A4^A(r)i 
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Note that a (r) in (lAlSp satisfies the equation a (t) cr (r) — ct^ (t) + A (r) cr'* (t) = 0. Once we have obtained cr (r), 
we have 

^l.^ (r) = / Aia^ (t) dr + Aq and ^^^ (r) = / A2(J^ (r) dr + A7, (A16) 

where Ag and A7 are reo/ constants. Integrating, we get 

, . 2AaA-\ 1 , , , _> 

^-(^) = -^m 7 ^^T + ^6, (A17) 

V-4 i'^) 1 + exp (2^4\/^(r)TJ 
and 

A*,(r) = -^^ / ' . N +-47. (A18) 

V^4(0 1 + cxp (2A4 v/;4Mr j 

Assuming the following boundary conditions cr (0) = cto > 0, ^^{too) — /i^, (tqc) = 0, we find that rro ~ ^4, 
Aq = At = 0. Finally, the geodesic trajectories become, 

2(70^1 1 , ^ 2(70^2 1 



^A (r) 1 + cxp {2aQyfMr)T^ ' ^ ' ^^l (r) 1 + exp {^gq^/AJ^t^ 



cxp ( (To V-4 ('')'' 

a{T;T) = 2(70 ^7 ^, (A19) 

1 + exp ( 2(70 v/-4(r 



with A(r) defined in \Kh\ . A\ and A2 real constants and ctq > 0. Note that a (r; r) e (0, + 00) while /i^. (t; r) and 
A<y (r; r) G (-00, +00). 
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